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Abstract. We describe a method for the accurate calculation of bound-state and 
resonance energies for one-dimensional potentials. We calculate the shape resonances 
for symmetric two-barrier potentials and compare them with those coming from the 
Siegert approximation, the complex scaling method and the box-stabilization method. 
A comparison of the Breit-Wigner profile and the transmission coefficient about its 
maximum illustrates that the agreement is better the sharper the resonance. 



1. Introduction 

In a recent paper Rapedius[T] showed how to calculate resonance positions and widths by 
means of the Siegert approximation (SA). He applied it to two exactly solvable models 
and also to the transmission through a double barrier. In the latter nontrivial example 
he compared the approximate SA resonances with the more accurate ones provided by 
complex scaling (CS). In this interesting pedagogical paper Rapedius showed that the 
SA is suitable for narrow resonances. He also described the difference between Siegert 
and transmission resonances. The former are complex eigenvalues of the Schrodinger 
equation whereas the latter are related to the maxima of the transmission coefficient. 

Somewhat earlier Dutt and Kar[2] had discussed scattering through smooth double 
barriers constructed by means of Gaussian functions. They compared the accurate 
transmission coefficient calculated numerically with analytical expressions derived by 
means of the WKB method. The agreement is remarkable for all values of the energy 
of the incident particle. 

The physical interest in double-barrier potentials arises from models for fission 
barrier [31 H] as well as in the study of chemical reaction thresholds [5] and simple 
molecular collisions [6]. In such cases sharp resonances are associated to states trapped 
between the two barriers [U [5] and exhibit Lorentzian or Breit-Wigner (BW) profile[6]. 
In addition to its remarkably accurate results [1] the WKB method is suitable for proving 
that the resonances are of BW type[8], [9]. Square wells and barriers have also proved 
suitable for a pedagogical discussion of bound states, virtual states and resonances |10j. 
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The purpose of this paper is to illustrate the relation between Siegert and 
transmission resonances in a somewhat more detailed way by means of the symmetrical 
two-barrier potential discussed by Rapedius[T]. The Schrodinger equation for this 
potential is not exactly solvable but one can calculate the Siegert resonances accurately 
by means of the Riccati-Pade method (RPM)[m [T2] . It is also possible to obtain 
the transmission coefficient as a function of the energy by means of the Wronskian 
method [ini [Tl] and thus to compare its resonance peaks with the BW expression for 
different barrier heights. This paper is expected to be a useful complement to those 
earlier pedagogical discussions on the subject [H Ej [10] and a suitable approach to the 
problem for courses of quantum mechanics at advanced undergraduate or graduate level. 
In section[2]we derive the SA in a way that differs from that followed by Rapedius[T]. 
Like this author we introduce the Siegert and transmission resonances that will be 
discussed and compared in later sections. In section |3] we outline the main ideas of the 
RPM for symmetric one-dimensional potentials. In section |4] we calculate the Siegert 
resonances by means of the RPM and the transmission coefficient by means of the 
Wronskian method. We can thus compare the peak of the transmission coefficient with 
the BW profile for several barrier heights. Finally, in section [5] we summarize the main 
results and draw conclusions. 

2. The Siegert approximation method 

RapediusP developed the SA from the time-dependent Schrodinger equation. Since 
this equation is not used in the calculation of neither the Siegert nor the transmission 
resonances we think that it may be fruitful to derive the SA entirely from the time- 
independent Schrodinger equation. 

In order to simplify and facilitate both the algebra and the numerical calculations 
in this paper we first convert the Schrodinger equation into a dimensionless eigenvalue 
equation. In this way one removes all the physical constants and reduces the number of 
model parameters to a minimum. Although we have already discussed this well known 



Resonances for symmetric two-barrier potentials 4 

procedure in earlier papers [T^ [T^ [T5| we think that it is worthwhile to insist on the 
advantages of working with proper dimensionless equations in physics. 

The time-independent Schrodinger equation for a particle of mass m that moves in 
one dimension (— oo < X < oo) under the effect of a potential V{X) is 

- ^nX) + V{X)i;iX) = E^{X) (1) 

where a prime indicates derivative with respect to the coordinate X. If we define the 
dimensionless coordinate x = X/L, where L is an appropriate length scale (or length 
unit), then we obtain the dimensionless eigenvalue equation 

(p"{x) + v{x)ip{x) = eip{x) 

I— , -. , . niL"^ , , niL'^E , , 

^{x) = VL^(Lx), v{x) = -^V{Lx), e = -j^ (2) 

The length unit L that renders both e and v{x) dimensionless is arbitrary and we can 
choose it in such a way that makes the Schrodinger equation as simple as possible as 
shown in the first example of section HI 

For the time being we just assume that the potential tends to zero 

lim v{x) = (3) 

|x|— >oo 

faster than |x|~^ and that the particle approaches the interaction region from the left. 
Therefore, the boundary conditions are 

lim ip(x) = (p-(x) = Ac''"'' + Be''-''^ 

X— >— oo 

lim ip(x) = ^+(x) = Ce^*^^ (4) 

where k = \/2e. 

It follows from the Schrodinger equation (l2l) and its complex conjugate that 

(v.V*-V''V)' = 2(e*-e)|v,p (5) 

If e is real then the Wronskian W{(f*, (p) = (p'(p* — (f'*(p is constant and the substitution 
of the asymptotic solutions (j4]) into W{ip*i,ip-) = W{(p^,(p+) leads to 

\A{e)\'-\B{e)\' = \C{e)\' (6) 
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or r(e) + R{e) = 1 where T(e) = |C(e)|VI^(e)P and R{e) = |5(e)|V|A(e)|2 are the 
transmission and reflection coefficients, respectively. 

Following Rapedius[T] we call transmission resonance to an energy value e = e-r 
for which T(e) exhibits a local maximum. In the particular case that the potential is 
symmetric about the origin v{~x) = v{x) the maximum has unit value T{eT) = 1[8]; 
consequently, R{eT) = 0, |A(er)P = |C(eT)P and |5(er)P = 0. The boundary 
conditions for such particular energy value are 

V2t±(^) = ikT^T±ix) (7) 

and 

1,^ |2 I /I |2 ir^ |2 \,„ |2 fQ\ 

where the subscript T indicates that e = e^. For concreteness we assume that the 
potential is symmetric from now on. If it were symmetric about another point x = Xq 7^ 
we would simply shift the coordinate origin from x = to x = Xq. 

A resonance (or Siegert) eigenstate is a solution to the Schrodinger equation that 
behaves asymptotically as a purely outgoing wave [71 E] (and references therein): 

lim ifsix) = ifs-i^) = Bse-'^^'' 
x— >— 00 

jimy,5(x) =^s+{x) = Csd''- (9) 

where the subscript S indicates that e = e^es which is a complex eigenvalue e^es = ^R+i^i- 
The real part is the resonance position and the imaginary part is related to the resonance 
width r = — 2e/ > 0. For practical purposes it is customary to assume that these 
boundary conditions are approximately valid for a sufficiently large coordinate value 
|x| = a. Therefore, if we integrate equation ([5]) between —a and a we obtain 

(V^V* - ^'*^)W = ^(ks + k*s) (IV.5+P + \^s-\') 

= -Aiei r \^s?dx (10) 

J —a 
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The Siegert states for a symmetrical potential are either even or odd; therefore 
|(/95+P = |(y95_p. In the case of a sufficiently narrow resonance |e/| <^ e^ we may 
safely carry out the additional approximation that e^es ~ ^r ~ ^t- 

Since the Siegert state is strongly localized in the well between the barriers located 
at a; = ±6 [1] we can also write 

/a rb rb 

\ips\^dx^ / \ips\^dx^ / \^T? dx (11) 

-a J~b J-b 

Finally, from the equations above we obtain an expression for the resonance width: 

r = ^^^M^ (12) 

ro I 19 J ^ '^ 

Jo l^rrdx 
already derived by RapediusP by means of the time-dependent Schrodinger equation. 

A more rigorous, general and elegant derivation of the results above was given by 
Whitton and Connor [T6] many years ago by means of a Wronskian analysis. Those 
authors called the S and T boundary conditions "outward moving waves only" and 
"forward moving waves only", respectively. Rapedius[T] also showed an expression 
for the nonsymmetric case v{—x) ^ v{x) that we do not consider here. The general 
expressions of Whitton and Connor [16] do in fact apply to both the symmetric 
and nonsymmetric case. Note that the energies for the forward-moving-waves-only 
boundary conditions are real for a symmetric potential but they may be complex for a 
nonsymmetric one [T6] . 

It is clear from all the assumptions made above that the SA equation ( IT2i) applies 
only to sufficiently narrow resonances. In fact, Rapedius [1] illustrated this point by 
comparing SA results with those provided by the more accurate CS. In section H] we 
propose an alternative comparative discussion of Siegert and transmission resonances. 

3. The Riccati-Pade method (RPM) 

In order to compare Siegert and transmission resonances we need sufficiently accurate 
complex eigenvalues of the Schrodinger equation with purely outgoing-wave boundary 
conditions. There are many suitable methods for this purpose; for example, Rapedius[l] 
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discussed the SA and CS ones. In what follows we outline the RPM that yields 
remarkably accurate results for narrow and broad resonances [TTl [T2] (a pedagogical 
approach to RPM for bound states is available elsewhere [17]). We think that the 
RPM may be a suitable practical tool for a quantum mechanics course because the 
derivation of its main equations and their implementation in a computer program are 
both straightforward. 

Suppose that we want to obtain the eigenvalues of the dimensionless Schrodinger 
equation ([2]) with a symmetric potential v{—x) = v{x) (without loss of generality we 
assume that v{0) = 0). We restrict ourselves to those eigenstates that are either even 
^p{—x) = (p{x) or odd (f{—x) = —^p{x) (both bound and resonance states satisfy this 
criterion) . We define the regularized logarithmic derivative of the eigenfunction 

where s = or s = 1 for even or odd states, respectively. It satisfies the Riccati equation 
fix) + -fix) - fix)' + 2t;(x) - 2e = (14) 

X 

If we can expand w(x) in a Taylor series about x = 

oo 

f(x) = ^fjX^-' (15) 

i=i 

then we can also expand fix) about the same point as 

oo 

fix) = xY.f,ie)z^, z = x' (16) 

i=o 

Note that the term s/x in equation ( fT3l) removes the pole at origin in the case of odd 
states iip°'^'^iO) = 0). The expansion of equation flT^ in a Taylor series about the origin 
leads to a recurrence relation for the coefficients fj that enables us to obtain as many 
coefficients /,(e) as necessary: 
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Since f{x) has poles at the zeros of (p{x) we look for a rational approximation of 
the form f[x) ~ x[M/N][z), where 

y-Af j M+N+1 



Because we can arbitrarily choose bo = 1 we are left with M + N + 1 coefficients aj 
and bj of the rational function and the unknown energy as independently adjustable 
parameters. Therefore we require that the rational approximation (Pade approximant) 
yields M + N + 2 exact coefficients of the Taylor series for f{x) as explicitly indicated 
in equation (TTSl) . When M > N we easily derive the following equations: 



ram{j,N) 

J2 fefc/i-fc = flj, j = 0, 1,...,M 

fc=0 

N 

Y, bkfj-k = 0, J = M + 1, M + 2, . . . , M + iV + 1 

fc=0 



(19) 



We can view the second set of equations as a system of A^ + 1 homogeneous equations 
with A^ + 1 unknowns bNibN-ii ■ ■ ■ ,bo. Therefore, there will be a nontrivial solution 
only if e is a root of 

fu-N+l fM-N+2 ■ ■ ■ Jm+I 
fM-N+2 Im-N+S ■ ■ ■ /m+2 



Hiie) 







(20) 



/m+I fM+2 ■ ■ ■ fu+N+l 

where d = M — N = 0,1,... and D = A^ + 1 = 2,3,... is the dimension of the 
Hankel determinant Hf)(e). As D increases, sequences of roots e^^''^ of the Hankel 
determinant converge toward the allowed energies of the Schrodinger equation. This 
sort of quantization condition applies to bound states and resonances |1H [T2| [T7] . We 
simply identify convergent sequences of real and complex roots e^^''^ of the Hankel 
determinants H'lj{e) for D = 2, 3, ... , D^ and estimate the error of the calculation as 



w 



[D„ 



w 



[Dm-1] 



where w'^' stands for the real or complex part of e'^''^'. We thus 



truncate the results to the last stable digit. 
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4. Results and discussion 

As a first illustrative example we consider the symmetrical double barrier V{X) = 
VoX'^e~°'^ , where Vb,a > 0, discussed by Rapedius[l]. The Schrodinger equation for 
this potential is not exactly solvable but it is sufficiently simple for pedagogical purposes. 
In this case the dimensionless potential is given by 

vo = 2 ^ A = aL (21) 

Note that without loss of generality we can treat this problem as a one-parameter model 



because we may have either vq = 1 when L^ = h/y/mVo or A = 1 when L^ = 1/a. 
However, we write v{x) as a two-parameter potential following Rapedius[T]. This 
potential exhibits a well centered at x = and two barriers of height vi, = vo/{eX) 
symmetrically located at x = ±6, b = 1/vA- 

For comparison purposes we first consider the potential parameters fo = 1/2 and 
A = 0.1 already chosen by Rapedius[I]. Table [U shows the first six resonances obtained 
from Hankel sequences with D < 20 and d = 0. The first three of them agree with those 
calculated by means of SA and CS[I]. Present results are supposed to be accurate to 
the last digit and, consequently, much more accurate than the SA and CS ones shown 
by Rapedius[T]. However, it is worth mentioning that the CS results can in principle be 
made as accurate as desired by increasing the dimension of the basis set [19]. As stated 
above (and already shown in Table [1]) the RPM is suitable for both sharp and broad 
resonances. However, the accuracy of the RPM for a given determinant dimension D 
decreases with the resonance "quantum number" n because the Hankel sequences for 
higher resonances appear at greater values of D. 

The RPM is extremely accurate but it is not as general as other approaches (like, 
for example, CS) because it only applies to separable Schrodinger equations. However, 
one-dimensional and separable models are widely discussed in most courses on quantum 
mechanics. 
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It is not difficult to calculate the transmission probability T(e) for the potential 
(I2T]) by means of the Wronskian method [T3l ITi] and thus compare the Siegert e^es and 
transmission e^ resonances discussed by Rapediusjl]. More precisely, we compare the 
actual transmission coefficient T(e) and the BW shape 

in a neighbourhood of the maximum of T(e) a.t e = ej- ^ e^i. Connor [S] has already 
shown that T{eT) = 1 for a symmetric potential and that equation (1221) is a reasonable 
approximation for isolated sharp resonances. 

Table [2] shows the lowest resonance for vq = 2, 5, 10, 15 and A = 1 calculated by 
means of the RPM exactly as discussed above. We have rounded off all the results to the 
ffist unstable digit. The second column displays the approximate values of the barrier 
heights Vb = Vq/c. Note that \vb — e^l increases with t>o; that is to say, the resonance 
moves deeper into the well between the barriers as Vq increases. At the same time it 
becomes narrower and therefore more stable (its decay rate F = — 2e/ decreases). Fig. [T] 
shows the potential and the location of the lowest resonance (horizontal line) for each 
case and illustrates graphically the behaviour just discussed. 

Fig.|2]shows the transmission probability T(e) calculated by means of the Wronskian 
method [131 Ej and the BW profile (1221) with the values of e^ and e/ given in Table |2j 
It clearly illustrates that T(e) becomes sharper about the maximum as F decreases. 
In addition to it, we appreciate the well known fact that the BW proffie gives a better 
description of the peak of the transmission coefficient the sharper the resonance. In other 
words, as the potential parameter increases from fg = 2 to fo = 15 and F decreases 
the BW profile (122!) becomes increasingly more accurate in the neighbourhood of the 
maximum. 

The BW profile is suitable for isolated resonances, and closely spaced broad 
resonances may overlap. In such a case T(e) becomes rather too distorted for the BW 
profile to fit satisfactorily. Fig. [3] shows the two first resonances for several values of wq. 
The resonance positions are indicated by symbols and their widths by error bars. We 
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clearly appreciate that the resonances overlap for f o = 2 and tig = 3 in which cases the 
BW profile is unsuitable. The magnitude of the overlap diminishes with vq as the gap 
between the resonances increases and their widths decrease. 

Klaiman and Moiseyev[T8] have recently proposed an improved profile based on the 
same information required for the BW one. The KM profile is nonsymmetric and it is 
peaked at le^esl > ^r which corrects the fact that typically eji < ex as shown in Fig. |2j 
We do not show the KM profile here because the correction is mild for the present model. 
The reason is that ej <^ e^ and le^esl ~ ^b. for the isolated resonances (say fo > 4). 

Once we have the complex energy eigenvalue we can easily calculate the Siegert 
state by numerical integration or any other approach. For example, its Taylor expansion 
about the origin 

oo 

^ix) = j:c,x'^^' (23) 

i=o 

provides a suitable approximate analytical expression. It is not difficult to verify that 
the coefficients of this series are given by the recurrence relation 



Cj+l 



k=0 



J = 0,1,... (24) 



(2j + s + l)(2j + s + 2) 
where we arbitrarily choose cq = 1 . The Taylor expansion fl23l) converges for all values of 
X but in practice we can only add a finite number M of terms and the resulting partial 
sum for |c/)(x)p tends to infinity as x — > oo. For this reason the partial sum will be valid 
only in a finite interval —xm < x < xm that we should choose judiciously. Fig. |4] shows 
v{x) for Vq = 1/2 and A = 0.1 and |y9(x)p for M = 24. This figure is similar to the one 
shown by Rapedius [1] for the same potential on a wider abscissas interval. We clearly 
see that the Siegert state is so strongly localized in the well between the two barriers 
that |(y9(x)p ~ for |x| > b which justifies one of the assumptions made in section [2l 

Some time ago Korsch and Gliick[T9] proposed a pedagogical approach to the 
calculation of the eigenvalues of the Schrodinger equation by means of the matrix 
representation of the coordinate and momentum. Among other illustrative examples 
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they considered the one-dimensional potential 

v{x) = (y - ^) e^'^' + J (25) 

which supports one bound state and many resonances when J = 0.8 and A = 0.1. Note 
that this potential-energy function reduces to a particular case of Eq. ( 1211) when J = 0. 
Those authors calculated the bound-state energy by a straightforward application of 
the matrix method and the first two resonances by means of the box-stabilization 
method[19]. The application of the RPM to this problem is straightforward and 
one simply looks for converging sequences of real and complex roots of the Hankel 
determinants in order to obtain the energies of the bound-state and resonances, 
respectively. Table [3] shows the eigenvalues estimated from Hankel determinants of 
dimension D < 20. Also in this case we see that the RPM results are remarkably 
accurate. 

5. Conclusions 

One of the purposes of this paper is to compare the Siegert and transmission resonances 
in an alternative way to that discussed by Rapedius[l]. To this end we calculated 
accurate Siegert resonances by means of the RPM and the transmission coefficient by 
means of the Wronskian method [T3l [T3] for the two-barrier potential (I2T!) . The chosen 
model is nontrivial but sufficiently simple for the straightforward aplication of both 
approaches. In this way we can compare the shape of T(e) about its maximum at 
e = ej- (transmission resonance) with the BW profile fl22l) constructed by means of the 
complex eigenvalues eres = ^r + i^i (Siegert resonance). The comparison illustrates 
the well known fact that the BW profile fits sharp isolated resonances more accurately. 
This conclusion complements Rapedius observation that the SA becomes less accurate 
as the resonance width increases[T]. In addition to it, we have argued that in the 
present example the KM profile [T8] mildly improves on the BW one because the lowest 
resonance is rather too sharp when it can be considered sufficiently isolated from the 
next resonance. 
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We have derived the SA from the time-independent Schrodinger equation in a 
way that is more consistent with the methods for the calculation of the resonances 
described in this and Rapedius' papers. The Wronskian approach is most suitable for 
the discussion of the problem and enables the decomposition of the resonance widths 
into the sum of partial widths[16] which we have omitted here. 

We have shown that the RPM is a powerful tool for the location of resonances in 
one-dimensional quantum-mechanical problems (more precisely, it applies to separable 
Schrodinger equations). The RPM is a local approach based on the expansion of the 
logarithmic derivative of the wavefunction about a chosen point, and for this reason 
it may not be the most convenient method for a pedagogical illustration of the main 
features of the Siegert states (for example, the discussion of the boundary conditions). 
However, it provides remarkably accurate results and the derivation of its main equations 
offers no difficulty. In addition to it, the recipe for its application is simple and easy to 
encode into a computer program, especially in most available computer algebra systems. 
We think that the RPM may be a valuable tool for teaching purposes because it allows 
the students to obtain accurate bound-state and resonance energies and compare them 
with those provided by other approaches like SA, CS or the matrix methodp, [T9] . 
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Table 1. Lowest resonances for the potential ([21]) with vq = 1/2 and A = 0.1 



n 



e/? 



-e/ 




1 
2 
3 

4 
5 



0.46014727653933356360 9.6203883198201929683 x 10"^ 

1.2804203534682821470 1.6737132594145830404x10-^ 

1.8531086351750533910 6.7240255103872613345 xlO-^ 

2.2323252762455511600 0.33989855689185650713 

2.567615869399468602 0.8194028131702960163 

2.887957554267041665 1.409344599863779927 



Table 2. Lowest resonance for the potential ((2T|) and several potential parameters 



Vo 



Vb 



(^R 



e/ 



2 0.7358 0.55937118458252732995 0.15830525114271135525 

5 1.839 1.1082157629920295074 0.078972583905329832058 

10 3.679 1.7816763825869113601 0.023794309337967155927 

15 5.518 2.3042519331774868362 0.007347829662205245864 







Figure 1. (Color online) Potential function (blue) and quasi-stable energy (red) for 
vo = 2, 5, 10, 15 
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Table 3. Bound state and resonances for the potential ([JF 



Sr 



-£l 



0.5020403621419 

1.4209709457146932076 

2.1271970775224959319 

2.5845828598531001914 

2.9244219292377372486 

3.255486140023381540 

3.5572161626513698 

3.824329026868890 

4.055433668209184 

4.249963938764321 

4.407748386304 

4.528814027868 





5.82652808855403 x 10"^ 

1.5447312841805183109x10-2 

0.17375071916219928095 

0.564794965582576499 

1.1115316000246994816 

1.7555062346769250 

2.4874451532278992 

3.29864420145319 

4.18316582758871 

5.136439406966 

6.15480966701 






1 2 


3 4 5 



Figure 2. (Color online) Numerical transmission probability (solid line, blue) and 
Lorentzian profile (dashed line, red) for vq = 2, 5, 10, 15 
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Figure 3. (Color online) Real (symbols) and imaginary (error bars) parts of the first 
(blue, circles) and second (red, squares) resonance 




Figure 4. (Color online) Potential (blue) and |v3(x)p (red) for vq — 1/2 and A — 0.1 



